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First  Interim  Reoorc  (OOOiA.1)  on  the  work  performed  by  Numerirai  Computation  Cor 
during  the  period  of  Decemoer  i5,  1935  -  February  lO,  1936. 
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The  scattering  of  normal  incident  time-harmonic  TEM  electromagnetic  wave 

by  a  cylindrical  target  with  axis  along  z-direction. is  considered,  e.g., 

E  =  E  i  +  E  i  and  ii  =  H  i  ,  where  i  is  the  unit  vector  in  the  a-direction. 
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The  whole  space  domain  Q  is  divided  into  three 
connected  but  non-overlapping  sub-domains, 
the  interior  region  .1,  representing  the 

t 

target  and  possessing  a  non-orthcgonal 
cvlindrical  grid  system  centered  in 
itself,  the  intermediate  region 
representing  the  free  space  just  outside 
of  the  target  and  possessing  the  same 
grid  svstem,  and  the  exterior  region 
representing  the  tar-field  free  space 
but  truncatec  at  a  large  distance  away 
from  the  target  and  possessing  the 
standard  orthogonal  cylindrical  grid 
system  (Fig.  1) . 

To  facilitate  the  discretization 
of  the  Maxwell's  equations  cn  the  non- 
orthogonal  grid  system,  the  following 
integral  forms  of  the  Maxwell's  equations  are  used 

dj.  =  i.Aj  •  d£, 

<t  H  ,•  dl  -  B(j  -  i'A)e)*E,  •  ds . 

<j  £_•  dl  =  iuj  C^'J,  H  •  ds , 

®  H-*  dl  =  -iujE  E,  •  ds , 
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ct-  E^*  dl  =  i'^u  •  ds , 
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H  •  dl  =  'tS  (,  J  -  iojE  E  J 
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and  the  asymptotic  terminating  condition  at 

n  X  E-  =  (y  /€  )'^('n  x  h  )  , 
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where  n  is  the  unit  normal  vector  at  the  interfaces  and 
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Eq.  (.1)  is  discretized  by  using  the  rectangle  rale  on  the  line  integral 
around  the  edges  of  ail  incremental  quadrilateral  defined  by  the  grid  system. 
Let  each  grid  point  of  the  non-orthogonai  polar  grid  system  be  denoted  by 
(r..,  3.)  s  (i,j),  where  "i"  and  "j"  denote  the  i-th  closed  cylindrical  grid 
line  and  the  j-ch  radial  grid  line  respectively;  let  the  center  of  the 
quadrilateral  defined  by  (i,j),  (i+l,j),  ^,i,j+l)  and  (i+l,j+l)  be  denoted  by 
(i-r^,j+h);  let  incremental  distance  between  the  points  :x  and  3, 

and  be  the  area  of  the  incremental  quadrilateral  centered  at  y. 

Moreover,  let  i  -  1,2,3,...,!,  and  j  =  1,2,3,. ..,J. 
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ir.al  discretized  (1)  is 


H  .  1  .  .  j  — A  .  j  ,  .  • 

i+h,J+h’ 
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In  this  way,  the  most  natural  finite  difference  discretization  of  the  Maxwell's 
equations  on  a  staggered  grid  system  (Fig.  2)  is  obtained. 

If  the  differences  of  the  material  properties  spread  linearly  across  a 
grid  zone  instead  cf  across  the  interface,  then  there  is  no  need  to  impose  the 
boundary  conditions  (2)  at  the  material  interface,  because  the  boundary  condi¬ 
tion  for  the  tangential  component  of  ^  is  satisfied  automatically  and  the  other 
three  boundary  conditions  are  also  satisfied  automatically  but  approximately. 

In  this  wav,  there  is  no  cumbersome  instruction  and  treatment  at  the  interface 
to  slow  down  Che  calculation  on  the  computei'.  The  discretization  of  the 

terminating  condition  (3)  is 

E.  =  (h  /£  )'^  Iv  1  .  M  , 

I,3'''4  o  o  1-4, 3+4 

3  =  1,2,3, . . . ,J. 


At  this  moment,  we  have  completely  discretizated  (l)-(3)  according  to  the 
above  description.  We  are  Crying  to  organize  these  complex  algebraic 
equations  for  the  purpose  of  programming. 
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